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DETECTING A SIGNIFICANT EVENT IN 
EXPERIMENTAL DATA AND USE OF SUCH FOR 
DETECTION OF ENGAGEMENT DURING A 
5 MECHANICAL TEST 

The present application is based on and claims the 
benefit of U.S. provisional patent application Serial 
No. 60/444, 446, filed February 3, 2003, the content of 
which is hereby incorporated by reference in its 
10 entirety. 

BACKGROUND OF THE INVENTION 
Sensitive, accurate, and consistent detection of 
significant events in experimental data is an 
important capability for a myriad of applications. 

15 For example, in experiments or tests that measure 
mechanical properties of materials, it is often 
important to measure a significant event known as the 
engagement point. For contact-type tests used 
typically to measure a mechanical property such as 

20 hardness, and which typically utilize an indenter, 
the engagement point is the point at which the 
indenter first contacts and begins to apply force to 
the test surface. For tensile tests measuring, for 
example, stress and strain, the engagement point is 

25 the point at which all slack has been removed from 
the test apparatus and test sample, and any further 
movement in the same direction begins to stress and 
strain the test sample. In another example, it can be 
desirable to detect changes in temperature from a 
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baseline in a particular process or environment. 
Often it is difficult to detect a significant event 
using a simple threshold, especially when measurement 
changes are gradual. A system and method of 
5 accurately, sensitively, and consistently detecting 
significant events in experimental data would have 
significant utility. 

SUMMARY OF THE INVENTION 

10 One aspect of the present invention includes a 

computer readable medium that includes computer 
executable instructions for detecting a significant 
event in measurement data. Measurements are processed 
to calculate values indicative of multiples of a 

15 standard deviation. The point of onset of the 
significant event is calculated as a function of the 
calculated values. 

Another aspect of the present inventions also 
includes detecting the significant event by receiving a 

20 data series indicative of test measurements as a 
function of a variable. At least one processed series 
is generated from the data series. A first point and 
a second point are identified on at least one of the 
processed series. The significant event is then 

25 calculated based on at least one of the two points. 

Other aspects of the present invention also 
includes associated methods and testing systems 
adapted to determine a significant event using 
methods of the present invention described herein. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
FIG. 1 is a schematic diagram of an environment 
where the present inventions can be practiced. 

FIG. 2 is a schematic diagram of a computing 
5 environment useful in practicing the present 
inventions. 

FIG. 3 illustrates a block diagram of an 
embodiment of the present inventions and illustrates 
broad concepts helpful in understanding the present 
10 inventions. 

FIG. 4 illustrates a block diagram of an aspect of 
the present invention. 

FIG. 5 is a flowchart illustrating a method in 
accordance with an aspect of the present invention. 
15 FIG. 6 is a graph of measured force as a function 

of displacement. 

FIG. 7 is a graph of a second derivative of the 
graph illustrated in FIG. 6. 

FIG. 8 is a graph of processed signals as a 
20 function of time generated in accordance with the 
present inventions . 

FIG. 9 is an enlargement of FIG. 6 and illustrates 
a point corresponding to a significant event calculated 
in accordance with the present inventions. 

25 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
Referring now to FIG. 1, there is shown a 
schematic diagram of a testing system, such as an 
indenter testing system used to measure mechanical 
30 properties of materials. Testing system 100 is 
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modified in accordance with the present inventions 
for testing a sample of material 105. 

Sample 105 is placed at a known location on 
computer controlled X-Y table 107 with the surface to 
5 be tested facing up. Electromagnetically driven 
indenter arrangement 109 is provided, which is 
positioned over the sample 105. Indenter 109 includes 
current driven load coil 111 activated by the 
application of electrical current from computer 

10 controlled variable current source 117 to move the 
probe tip 115 downward into contact with sample 105. 
Once tip 115 contacts sample 105, a pre-selected force 
pattern is applied to indenter 109 by the programmed 
variation of the current applied to drive coil 111. In 

15 some embodiments, probe tip 115 can be in the form of a 
typical triangular pyramidal diamond probe with an end 
radius of about 500 Angstroms. 

Current source 117 is controlled by system 
computer 119, which also controls X-Y table 107. 

20 Displacement of probe 113 is measured by capacitive 
displacement gage 121, whose output is connected to DC 
displacement detector 123. Detector 123 digitizes the 
DC displacement signal, which is fed through a first 
digital voltmeter 125 to an input of computer 119 for 

25 further processing in accordance with the present 
inventions. Voltmeter 125 can provide a calibrated 
readout of probe displacement to an operator during 
testing procedures, if desired. 

Force applied to sample 105 through indenter probe 

30 115 is monitored by DC current detector 127, which 
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senses the DC drive current applied to load coil 111. 
The DC load current is digitized by detector 127 and 
fed through a second digital voltmeter 129 to a further 
input of computer 119 for further processing according 
5 to the present inventions. In some embodiments / 
computer 119 can be connected to a mass storage device 
131 where data and system operating parameters are 
stored. 

Using testing system 100 as described above, a 

10 sample 105 is positioned at a known location on X-Y 
table 107 and programmed computer 119 is signaled to 
start a timed test and data gathering procedure. 
Generally, computer 119 is programmed to perform one 
or more prescribed indentation tests by automatically 

15 causing single or multiple indentations at designated 
locations on sample 105. Probe 115 is lowered at a 
very slow rate until contact is made with sample 105. 

During the test procedure, computer 119 records 
measurement data, such as force taken from DC current 

20 (load) detector 127 and probe displacement taken from 
DC displacement detector 123. Measurement data can be 
used later to determine mechanical properties of sample 
105. Additionally, measurement data such as values 
associated with measured force, displacement, and time 

25 can be processed using the present inventions to 
determine significant events occurring during the test 
procedure, such as when probe tip 115 "engages" sample 
105 during the indentation procedure. 

FIG. 2 and the related discussion provide a brief, 

30 general description of a suitable computing environment 
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in which the present inventions can be implemented. 
Although not required, testing system 100 and the 
present inventions can be described, at least in part, 
in the general context of computer-executable 
5 instructions, such as program modules, being executed 
by computer 119, 219. Generally, program modules 
include routine programs, objects, components, data 
structures, etc., which perform particular tasks or 
steps or implement particular abstract data types of 

10 the present inventions. Program modules and associated 
methods are described below and illustrated using block 
diagrams and flowcharts. Those skilled in the art can 
implement the block diagrams and flowcharts to 
computer-executable instructions. Moreover, those 

15 skilled in the art will appreciate that the inventions 
may be practiced with other computer system 
configurations, including multi-processor systems, 
networked personal computers, mini computers, main 
frame computers, and the like. The inventions may also 

20 be practiced in distributed computing environments 
where tasks are performed by remote processing devices 
that are linked through a communications network. In a 
distributed computer environment, program modules may 
be located in both local and remote memory storage 

25 devices. 

Computer 219 illustrated on FIG. 2 (similar to 
computer 119 illustrated on FIG. 1) , comprises a 
conventional personal or desktop computer having 
central processing unit (CPU) 232, memory 234 and 
30 system bus 236, which couples various system 
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components, including memory 234 to CPU 232. System bus 
236 may be any of several types of bus structures 
including a memory bus or a memory controller , a 
peripheral bus, and a local bus using any of a variety 
5 of bus architectures. Memory 234 includes read only 
memory (ROM) and random access memory (RAM) . A basic 
input /output (BIOS) containing the basic routine that 
helps to transfer information between elements within 
computer 230 , such as during start-up, is stored in 

10 ROM. Storage devices 238, such as a hard disk, a floppy 
disk drive, an optical disk drive, etc., are coupled to 
system bus 236 and are used for storage of programs and 
data. It should be appreciated by those skilled in the 
art that other types of computer readable media that 

15 are accessible by a computer, such as magnetic 
cassettes, flash memory cards, digital video disks, 
random access memories, read only memories, and the 
like, may also be used as storage devices. Commonly, 
programs are loaded into memory 234 from at least one 

20 of storage devices 238 with or without accompanying 
data. 

Input device 24 0 such as a keyboard, pointing 
device (mouse) , or the like, allows a user to provide 
commands to computer 230. Monitor 242 or other type of 
25 output device is further connected to system bus 236 
via a suitable interface and provides feedback to a 
user. 

Control signals 217 output to testing system 100 
comprising current source 117 and X-Y table 107 
30 (illustrated in FIG. 1) based on program modules 
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executed by computer 119, 219 and through suitable 
interface 244. Interface 244 receives output signals 
221 indicative of displacement from voltmeter 125 and 
signals 221 indicative of force from voltmeter 129. 
5 Response signals 221 are processed based on program 
modules executed by computer 119, 219 to calculate or 
approximate, for example, material properties or 
significant events occurring during the testing 
procedure, such as engagement point. 

10 FIG. 3 illustrates a block diagram of an 

embodiment of the present invention, which can be 
implemented with a computer such as computer 119, 219 
executing instructions stored on a computer-readable 
medium. Modules, sub-modules, and associated steps 

15 for carrying out such instructions are described 
below. Those skilled in the art will readily 
understand that the illustrated modules, sub-modules, 
and steps can be combined or divided as desired. It 
is noted, that all steps of the present "inventions 

20 can be performed after the test procedure or 
experiment is completed. Optionally, some steps can 
be performed during the actual test procedure. Also, 
variables x, y, z, t, etc. all represent series of 
values that are controlled, sensed or measured during 

25 the test procedure or experiment, or series that are 
calculated or derived from such. Examples of such 
series include force y, displacement x, time t, 
reading number, stress, strain, flow rate, 
temperature, and the like. It should be noted that x, 

30 y, z, and t may actually represent any experimental 
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series. It should further be noted that x, y, z, and 
t can each be independent or dependent series. For 
example, in a testing system that measures mechanical 
properties of materials, where force is controlled 
5 and displacement is a sensed response, y can 
represent force and t can represent displacement. 
Finally, it should be noted that each time the order 
(n) changes, the series represented by t may also 
change. 

10 In FIG. 3, broad modules or steps that are 

helpful in understanding the present inventions are 
illustrated. Testing system 100 (illustrated on FIG. 
1) generates signals indicative of measurements, 
variables, or data 302 such as x, y, z, t described 

15 above. Measurements 302 are received as input to 
significant event determination module 300, which 
comprises pre-processing module 304, determination 
module 306, and optionally, uncertainty module 308. 
Pre-processing module 304 receives and processes 

20 measurements 302 in accordance with the present 
inventions to generate values used by determination 
module 306 to determine or detect events that are 
significant, such as engagement point S of probe tip 
115 and sample 105 (illustrated on FIG. 1) . Notably, 

25 such an engagement point has values associated with 
variables such as force, displacement, and time at 
engagement. Significant event determination module 
300 generates measurements or values for the 
variables or parameters at engagement point S, and 

30 optionally, an uncertainty as indicated at 310. 
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FIG. 4 illustrates an embodiment of the present 
invention and generally illustrates greater detail 
than the embodiment illustrated in FIG. 3. For 
convenience, the last digits of the reference 
5 numerals in FIGS. 3 and 4 correspond where possible 
(such as 302, 402) to indicate similar items or 
steps. FIG. 5 illustrates a flowchart of an algorithm 
or method 500 of the present invention. The blocks in 
FIG. 5 are labeled with letters "A" through "P" 
10 generally indicating steps of one or more methods. 
However, it should be noted that the indicated steps 
and their respective order are illustrative only and 
can be combined, recombined, divided, or reordered as 
desired. 

15 Testing system 100 (illustrated on FIG. 1) 

generates values associated with or indicative of 
measurements or data 302, 402. Measurements 302 , 402 
are provided to or received by significant event 
determination module 400. FIG. 5 generally 

20 illustrates a step of receiving measurements as start 
block 502 and block 504 where the series order, n is 
initialized. 

Significant event determination module 400 
includes non-zero trend detector 422, which is 

25 adapted to or capable of detecting trends in 
measurements 302, 402. FIG. 5 illustrates a step of 
detecting a non-zero trend in received measurements 
at block 506. Significant event determination module 
400 can also include pre-processing module 404, shown 

30 in dotted lines, that corresponds to pre-processing 
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module 304 (illustrated in FIG. 3) . Detector 422 can 
use any known means such as mathematical methods, 
electronic methods, curve fitting, etc. to detect a 
slope, change, or non-zero trend in measurements 302, 
5 402, especially prior to the significant event to be 
detected. Measurements 302, 402 can include variables 
such as force recorded or measured as a function of 
at least one other variable such as displacement 
and/or time. In some embodiments, force is a 
10 controlled variable and displacement is a response 
variable being measured. 

Generally, non-zero trends can be caused by the 
physical environment or the testing system itself. 
For example, a load cell may have some known offset, 
15 for example, 5mN. Another example is a linear 
relationship between force and displacement, which 
can exist due to internal stiffness of testing system 
100. In the examples provided, both the offset in 
force of 5mN and the internal stiffness cause a non- 
20 zero trend before the onset of the significant event 
or engagement point. In the present inventions, an 
initial non-zero trend in a data measurement signal 
is processed by "zeroing" or generating at least one 
processed signal or series that scatters or 
25 fluctuates about zero. 

In contrast, if non-zero trend detector 422 does 
not detect a non-zero trend, i.e. measurements 302, 
402 are scattered about zero as a function of another 
variable, such as displacement or time, then 
30 measurements 302, 402 can be provided to integration 
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module 314 directly. However, if non-zero trend 
detector 422 detects a pre-signif icant event non-zero 
trend in measurements 302, 402, then measurements 
302, 402 are provided to zeroing module 424 for 
5 further processing. FIG. 5 illustrates a step of 
zeroing at block 508, which is indicated with dotted 
lines. The step of zeroing can be iterative as 
indicated at 510. The series order is advanced with 
each loop, pass, or iteration through zeroing module 

10 508 as indicated at block 512. 

Zeroing module 424 comprises compensator 410 
and/or differentiator 412 that are used to process or 
generate a new series or signal output that 
approximates or scatters about zero. Compensator 410 

15 or differentiator 412 is user selectable based on 
design criteria as illustrated at compensation/ 
differentiation selection block 514 illustrated in 
FIG. 5. In embodiments utilizing compensation, 
compensator 410, (assuming variable y is a function 

20 of t) constructs a new series by "compensating" 
values (generally, but not limited to, all values) of 
measurements 302, 402 using an analytic function 
¥(t) . FIG. 5 illustrates compensation at compensation 
block 516. The function W(t) can be known from the 

25 physical environment or it may be determined by 
methods such as curve fitting within some significant 
pre-event period using, for example, a low-order 
polynomial function or some other appropriate form. 

Therefore, in the above example where a load 

30 cell contains a 5mN offset or l P(t)=5mN, compensator 
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410 constructs a new or processed series such as yi(t) 
where yi(t) =y 0 (t) -V (t ) =y 0 -5mN . In the second example 
where a linear relationship exists between yo(t) and 
displacement t, then a linear fit to pre-event data 
5 produces the linear function ^(t) where ^(tHkt+b, k 
is slope, and b is the constant at t=0. Compensator 
410 constructs or generates a new series yi(t) where 
yi(t)«y 0 (t)-Y(t), or yi (t ) =y 0 ( t) -kt-b. 

In embodiments illustrative of differentiation, 

10 differentiator 412 takes successive derivatives, 
usually up to a first or second derivative to 
construct at least one new or processed series 
scattering about zero. Thus, in the above example 
where a load cell has a 5mN offset, series yo(t) is 

15 approximately equal to 5mN, and differentiator 412 
generates a new series yi(t) from y 0 (t) by taking a 
first derivative of y 0 (t) so that yi (t ) =y 0 ' (t ) =0 . 
Thus, the new series yi(t) now scatters about or is 
approximately zero. In the second example where a 

20 linear relationship exists between yo(t) and 
displacement t, then series yo(t) can be expressed as 
y 0 (t)=kt+b, where k is slope and b is a constant as 
above. On a first iteration through differentiator 
412 (illustrated as differentiation block 518 on FIG. 

25 5), differentiator 412 takes a first derivative of 
y 0 (t) to generate a new or first processed series 
yi(t), where yi (t) =y 0 ' (t) =k. This means that series 
yi(t) scatters about k after a first pass through 
differentiator 412. Differentiator 412 then loops and 
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generates a second derivative of yo(t), or a second 
processed series y 2 (t) where y 2 (t ) =yi' (t ) =y 0 ' ' (t) =0 . 

Thus, in the second example two loops or passes 
through differentiator 412 are needed to generate or 
5 construct series y 2 (t) scattering or fluctuating about 
zero. It should be noted that differentiation can be 
accomplished by any known means, including but not 
limited to numerical methods, electronic methods, or 
curve fitting. In many instances, it can be 

10 advantageous to use differentiation versus compensation. 
For example, with differentiation, it is not generally 
necessary to know the underlying functionality of 
series variable y n , e.g. it is not necessary to know 
the load offset, or precise linear relationship 

15 between y and t. It should also be noted that there 
can be a pre-event quadratic, cubic, or higher order 
polynomial relationship between y and t, which would 
seemingly call for third, fourth, and higher order 
derivative to generate a processed series scattering 

20 about zero. However, taking successive derivatives 
generally results in greater signal instability, which 
effectively limit the number of derivatives that can 
be taken. 

In another embodiment, both differentiation and 
25 compensation can be used to zero input measurements. 
Thus, assuming there is a quadratic relationship 
between two variables such as force y and displacement 
t, then initially, y 0 can be expressed as yo=at 2 +bt+c 
where a, b, and c are constants. Taking iterative 
30 first and second derivatives of yo yields 
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Y2 (t) =y' ' (t) =c- Compensation can then be used where 
W(t)=c so that new series y3(t)=0. Thus, compensation 
and differentiation can be combined, especially to 
limit signal instability when a derivative greater 
5 than a second derivative is needed to zero input 
measurements . 

Zeroing module 424 outputs or generates processed 
series y n also called zeroed or "derivative signal" 
413 where at least a portion of the derivative signal 

10 scatters about zero, especially pre-event data. It is 
important to reiterate, however, that there are 
situations where input measurements 302, 402, 
especially pre-event data, scatter about zero and are, 
therefore, not in need of compensation or 

15 differentiation. Thus, derivative signal 413 is 
intended to include a series, new or original, that 
scatters about zero with or without differentiation or 
compensation. Derivative signal 413 is provided to or 
received by optional integration module 414 for 

20 further processing. The flowchart of FIG. 5 
illustrates that integration is optional at block 520. 
However, in most situations integration is needed to 
implement the present inventions. In embodiments where 
integration is selected, integration module 414 

25 generates or constructs new series y n +i(t) also called 
"integral signal" 415 by integrating series y n with 
respect to t over a selected increment or period r . 
FIG. 5 illustrates the step of integration at 
integration block 522, which advances the series order 
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at block 524. The processed series y n +i(t) generated by 
integration module 414 can be computed as follows: 

J>Vi(0 = jy n (*)dt. Eq. 1 

t-T 

Thus, the series y n +i(t) generates values indicated at 
5 416 for any value of t, It is noted that integration 
can significantly increase sensitivity to significant 
event onset, especially over some previous period r, 
because gradual changes are cumulative or accumulated. 

In integration module 414 a period r can be 

10 selected so that r is several times greater than the 
period in which detection is desired. For example, if 
t represents time and detection is desired within 1 
second, then r could be selected around 5 seconds. If 
r is set too large, the small offsets in y n may be 

15 accumulated, and y n +i(t) may not be centered about 
zero . 

Optionally, integration module 414 can compute 
values 416 using the following equation: 



20 



^ + i(') = - \y n (0dt 9 Eq. 2 



which can be viewed as the mathematical expression of 
mean values of y n (t) over the previous period, r. 
Integration can be accomplished by several known 
means, including but not limited to, numerical 
25 methods, electronic methods, or curve fitting. 

Integral signal 415 having values 416 generated 
by integration module 414 are then output to or 
received by normalization module 418, which is 
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indicated at normalization block 526 on FIG. 5. It 
should be noted that normalization module 418 
generates the same results on output signal 417 
regardless of whether above Eq. 1 or Eq. 2 is used in 
5 integration module 414. Normalization module 418 
comprises standard deviation computation module 420, 
which computes or generates standard deviation a 421 
for series or values 416 of integral function 415 over 
some period (typically a significant interval) before 

10 the significant event. FIG. 5 illustrates computing 
standard deviation a 421 at block 528. 

Normalization module 418 processes or generates a 
new series y n +i(t) by dividing all values 416 in series 
y n by standard deviation a as indicated at 422. FIG. 5 

15 illustrates the step of dividing y n by the standard 
deviation a at block 530 to advance the series order as 
indicated at block 532. Thus, the output of 
normalization module 418 are values 424 that are 
dimensionless and multiples of the standard deviation 

20 a . These values 424 can then be used to detect events 
such as the point of engagement. Thus, significant 
events are calculated or determined as a function of 
these values 424 that are indicative of standard 
deviation, especially multiples of the standard 

25 deviation. 

Normalization module 418 outputs values 424 that 
are received by determination module 306, 406. FIG. 5 
illustrates the broad step of event determination at 
block 534 indicated as dotted lines and comprising 

30 several sub-blocks 536, 538, and 540. Engagement 
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determination module 306, 406 comprises user selected 
threshold settings hi and h 2 as indicated at 426. The 
threshold setting hi is selected as a lower bound to 
indicate that the significant event has not yet 
5 occurred. The threshold setting h 2 is selected to 
indicate an upper bound or a point where the 
significant event is assumed to have already 
occurred. It is generally desirable to set threshold 
setting h 2 in consideration of a desired confidence 

10 level. Generally, if the pre-event distribution of 
y n (t) is approximately Gaussian, then a threshold 
setting h 2 in the range of 10 and 20 (i.e. where 
values 424 are between 10 and 20) yields good 
results. As with h 2r threshold setting hi should be 

15 set in consideration of the desired confidence level. 
Generally, selecting a threshold setting hi at a value 
less than 1 (i.e. where values 424 are less than 1) 
yields good results. 

Engagement determination module 406 then uses 

20 threshold settings h 2 and hi to determine points S 2 
and Si respectively, as indicated at 428. The point S 2 
is determined using a processed series y n (t) and 
threshold h 2 . S 2 is the first point in the series 
y n (t) for which there is excellent evidence for the 

25 onset of the event. Thus, S 2 is the first point, 
actual or interpolated, for which y n (t) differs from 
zero by a threshold, h 2 . Thus, if the threshold h 2 is 
set at 15, then the point S 2 occurs where the value of 
y n (t) first exceeds 15. Once S 2 has been determined, 
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it can be associated or correlated with points in 
other data series of the experiment. 

Point Si is then determined starting with the 
point S 2 and stepping back through values in the 
5 series y n (t) . Point Si is the first encountered point, 
actual or interpolated, for which y n (t) is within a 
given absolute-value threshold setting, hi. Absolute 
values are used because values for series y n (t) can be 
positive and negative in some situations. FIG. 5 

10 illustrates steps of determining points S2 and Si at 
block 536 and block 538, respectively. 

Determination module 306, 406 calculates an 
event such as engagement, herein denoted as point S 
or S' (illustrated in FIGS. 8 and 9) using S Xr S 2 , and 

15 any combination of data series. FIG. 5 illustrates 
the step of calculating S, S' at block 540. 
Calculation of S, S' may be accomplished in any 
number of ways. For example, it can be sufficient to 
set S to be the same point as S 2 , or S can be 

20 calculated to be the first actual or interpolated 
point after Si, for which y n (t) differs from zero by a 
less stringent threshold h s of, say, 4. It should be 
noted that in alternate embodiments, engagement point 
denoted as S' does not necessarily have to lie on a 

25 particular experimental trace z(t). For example, if z 
represents force and t represents displacement, it 
may be useful to calculate S' as the unique 
intersection point of linear extrapolations of z(t) 
data before Si and after S 2 (as illustrated as 

30 intersection point S' on FIG. 9) in which case, S' 
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may have associated values of force and displacement 
that are not actual measured values. Once S (or S' ) 
has been determined, it can be corresponded with 
values in other data series for the experiment, such 
5 as temperature. FIG. 5 illustrates the step of 
determining S at block 540. 

Significant event determination module 400 can 
include uncertainty module 308, 408 to calculate 
uncertainty Sz . The step of determining uncertainty is 

10 illustrated as block 542 on FIG. 5. Uncertainty Sz is 
the uncertainty in any sensed, response, controlled, 
or calculated series due to the uncertainty in the 
point of onset of the event. Calculation of Sz can be 
accomplished any number of ways. A conservative 

15 estimate of Sz might be obtained by calculating half 
the variation in z between Si and S2. Thus, 
uncertainty can be calculated as a function of Si and 
S 2 . If extrapolations are used to obtain S, 
uncertainties can be based on uncertainties in the 

20 fitted coefficients of the extrapolations. Significant 
event determination module 400 generates engagement 
point S and uncertainty Sz as indicated at 410, and 
illustrated at end block 544 on FIG. 5, that are then 
available to a user and/or for further processing. 

25 An example 

To demonstrate the present inventions, a hardness 
test was performed using a testing system similar to 
testing system 100 illustrated in FIG. 1. Although not 
required, such a testing system can include a computer- 

30 controlled AC signal generator fed to current source 



-21- 

117 and an AC displacement detector in line between 
capacitive displacement gage 121 and computer 119. Such 
a testing system has been referred to as a Nanoindenter 
XP, which is manufactured by MTS Systems Corporation of 
5 Eden Prairie, Minnesota and described in U.S. Patent 
No. 4,848,141, which is herein incorporated by 
reference. Although not illustrated, a similar test 
procedure can be performed on a tensile testing system 
such as described in U.S. Patent No. 6,679,124, herein 

10 incorporated by reference in its entirety. In the 
context of a hardness test, the "significant event" of 
interest is engagement between an indenter or probe tip 
115 and the surface of a test sample such as sample 
105. Force (y 0 ) , displacement (x) , and time (t) were 

15 recorded during the initial part of a hardness test on 
polycarbonate before, during, and after engagement or 
contact. Polycarbonate was chosen because it is a 
rather compliant material, which makes engagement 
detection more challenging. 

20 FIGS. 6-9 are pictorial graphs that illustrate 

results of a test procedure performed in accordance 
with the present invention. However, it should be noted 
that although pictorial graphs are helpful in 
understanding aspects of the present invention and can 

25 be rendered if desired to the user, pictorial graphs 
are not necessary for processing data associated with a 
significant event. FIG. 6 is a graph 600 that 
illustrates measurements of force (mN) as a function of 
a variable displacement (nm) . Both force, displacement, 

30 and time measurements were recorded during the test 
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procedure and received by a computer such as computer 
119, 219 (illustrated on FIGS. 1 and 2, respectively) 
for further processing in accordance with the present 
inventions. Graph 600 illustrates point S, indicated at 
5 602, which is the calculated engagement point according 
to the present inventions. To the left of S as 
indicated at 604, there is no engagement. To the right 
of S as indicated at 606, the probe tip is engaged with 
the sample. 

10 In order to obtain a series that is scattered 

about zero before contact, a plurality of new or 
processed series were generated by taking a first 
derivative and then a second derivative of force with 
respect to displacement. The second derivative was 

15 calculated using the following equation: 

ax 

It is noted that the notation y 2 (x) =y 0 ' ' (x) used above 
20 is equivalent to the notation in Eq. 3. FIG. 7 is a 
graph 700 illustrating the results for new series 
y 2 (x) plotted as a function of displacement x using 
Eq. 3. The differentiation was performed numerically 
using any point of interest and several data points 
25 backward (to the left) . In this context, repeated 
differentiation is advantageous because it can be done 
in real-time with no a priori information regarding 
the pre-engagement functional relationship between 
force and displacement. (In other words, one doesn f t 
30 need to know the offset or the stiffness of the 



-23- 

testing instrument.) The disadvantage is that repeated 
numerical differentiation of experimental data tends 
to produce series that are progressively more 
"scattered. " 

5 FIG. 8 is graph 800 illustrating series y 3 and y 4 

as a function of another variable time t rather than 
displacement x. In order to illustrate the advantage 
of integration such as illustrated by integration 
module 414 on FIG. 4 and block 522 in FIG. 5, 
10 processing was performed with and without integration. 
Thus, y 3 indicated at 802 is the series without 
integration and y 4 indicated at 804 is the series with 
integration. The series y 3 (t) at 802 was calculated 
with the following equation: 



15 



where a[y 2 (t)] is the standard deviation in y 2 (t) over a 
significant period before contact. 
20 In contrast, the series y 4 (t) is calculated by 

first computing a series, y 3 (t) (not shown on FIG. 8) 
with the following equation: x 



*(0= \y 2 (0dt. Eq. 5 

t-T 

25 

The integration was performed numerically using point- 
to-point trapezoid rule. The increment or period, r, 
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was set to 10 seconds. Using the y 3 (t) series 
generated from Eq. 5, the series y 4 (t) as indicated at 
804 is generated with the following equation: 



where cr[y 3 (t)] is the standard deviation in y 3 (t) from 

Eq. 5 over a significant time before contact. The 
series y 3 and y 4 indicated at 802 and 804 , 

10 respectively, are dimensionless . The values of y 4 
represent multiples of the standard deviation a from 
y 3 or integral values. It is apparent from FIG. 8 that 
the illustrated series y 4 more clearly indicates 
engagement point than series y 3 . For this reason, 

15 integration is generally desirable due to cumulative 
effects . 

Using the series y 4 (t) indicated as 804, the 
"upper" limit or bound S2 of engagement was determined 
to be the first point for which y 4 (t) exceeded 20. 

20 Thus, the threshold h 2 described above was set at 20 
(but could be set as low as 2 as desired) . The "lower" 
limit or bound for the point of engagement Si was 
determined by starting at S 2 and working backward to 
the first point where the absolute value of y 4 (t) is 

25 <1, i.e. threshold hi=l. The point S was acquired from 
the first point after Si where y 4 exceeded 4, i.e. 
threshold h s =4. All three points, Si, S 2 , and S are 
illustrated on FIG 8 as 808, 806, and 810 
respectively. FIG. 9 generally illustrates the data 
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illustrated on FIG. 6, but graph 900 is an enlargement 
of data or measurements around the point of engagement 
S indicated as 602. It is noted that points Si and S 2 
can be correlated across various series. For example, 
5 Si and S 2 correlate with 808, 908 and 806, 906 on FIGS. 
8 and 9, respectively. Engagement point S correlates 
with 602, 810 on FIGS. 6, 8 and 9. Finally, it is 
noted that point Si for which there is generally no 
discernable evidence of engagement or event onset can 

10 be obtained by any known means, including but not 
limited to the means described herein. The point Si, 
however determined, is then used in the calculation of 
engagement point S and/or uncertainty. 

Uncertainties corresponding to uncertainty 

15 module 408 and block 542 on FIGS. 8 and 9, 
respectively, were estimated for both force and 
displacement due to the uncertainty in the 
determination of the point of contact. A rather 
conservative estimate of uncertainty can be obtained 

20 by taking half the variation in both force and 
displacement between Si and S 2 . Thus, uncertainty can 
be calculated as a function of Si and S 2 . This method 
yielded uncertainty estimates of 0 . 6yN and 5.6nm, 
respectively. A more aggressive estimate of 

25 uncertainty is obtained by taking half the variation 
in force and displacement between the reading 
immediately before and after the point determined to 
be the point of engagement S. Thus, uncertainty can 
also be calculated as a function of S. This method 
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yielded uncertainty estimates in force and 
displacement of 0.077mN and 0.68nm, respectively. 

The present inventions, including the method 
presented in FIG. 5, are advantageous over simple 
5 prior art threshold detection because significant 
event onset (e.g. engagement point) in experimental 
data can be detected by more sensitive means. The 
detection means or apparatus also includes a means 
for determining lower and upper bounds of the range 

10 of possible event onset, or Si and S 2 . The present 
inventions also include a means for determining the 
uncertainty in any test parameter due to the 
uncertainty in the determination of S. Finally, prior 
knowledge of the pre-event functionality of the data 

15 series of interest is not generally required. 

Although the present invention has been described 
with reference to preferred embodiments, workers 
skilled in the art will recognize that changes may be 
made in form and detail without departing from the 

20 spirit and scope of the invention. For instance, 
although digitized signals are described herein for 
processing, analog signals and analog processing can be 
preferred as appreciated by those skilled in the art. 



